Measuring the Ellipticity of M 87* Images

The Event Horizon Telescope (EHT) images of the supermassive black hole at the center of the galaxy M 87 provided the first image of the accretion environment on horizon scales. General relativity predicts that the image of the shadow should be nearly circular, given the inclination angle of the black hole M 87*. A robust detection of ellipticity in the image reconstructions of M 87* could signal new gravitational physics on horizon scales. Here we analyze whether the imaging parameters used in EHT analyses are sensitive to ring ellipticity and measure the constraints on the ellipticity of M 87*. We find that the top set is unable to recover ellipticity. Even for simple geometric models, the true ellipticity is biased low, preferring circular rings. Therefore, to place a constraint on the ellipticity of M 87*, we measure the ellipticity of 550 simulated data sets of GRMHD simulations. We find that images with intrinsic axis ratios of 2:1 are consistent with the ellipticity seen from the EHT image reconstructions.Comment: accepted for publication to Ap

Axial focusing of impact energy in the Earth's interior: Proof-of-principle tests of a new hypothesis

A causal link between major impact events and global processes would probably require a significant change in the thermal state of the Earth's interior, presumably brought about by coupling of impact energy. One possible mechanism for such energy coupling from the surface to the deep interior would be through focusing due to axial symmetry. Antipodal focusing of surface and body waves from earthquakes is a well-known phenomenon which has previously been exploited by seismologists in studies of the Earth's deep interior. Antipodal focusing from impacts on the Moon, Mercury, and icy satellites has also been invoked by planetary scientists to explain unusual surface features opposite some of the large impact structures on these bodies. For example, 'disrupted' terrains have been observed antipodal to the Caloris impact basis on Mercury and Imbrium Basin on the Moon. Very recently there have been speculations that antipodal focusing of impact energy within the mantle may lead to flood basalt and hotspot activity, but there has not yet been an attempt at a rigorous model. A new hypothesis was proposed and preliminary proof-of-principle tests for the coupling of energy from major impacts to the mantle by axial focusing of seismic waves was performed. Because of the axial symmetry of the explosive source, the phases and amplitudes are dependent only on ray parameter (or takeoff angle) and are independent of azimuthal angle. For a symmetric and homogeneous Earth, all the seismic energy radiated by the impact at a given takeoff angle will be refocused (minus attenuation) on the axis of symmetry, regardless of the number of reflections and refractions it has experienced. Mantle material near the axis of symmetry will experience more strain cycles with much greater amplitude than elsewhere and will therefore experience more irreversible heating. The situation is very different than for a giant earthquake, which in addition to having less energy, has an asymmetric focal mechanism and a larger area. Two independent proof-of-principle approaches were used. The first makes use of seismic simulations, which are being performed with a realistic Earth model to determine the degree of focusing along the axis and to estimate the volume of material, if any, that experiences significant irreversible heating. The second involves two-dimensional hydrodynamic code simulations to determine the stress history, internal energy, and temperature rise as a function of radius along the axis

Evaluation of database technologies for the CTBT Knowledge Base prototype

This document examines a number of different software technologies in the rapidly changing field of database management systems, evaluates these systems in light of the expected needs of the Comprehensive Test Ban Treaty (CTBT) Knowledge Base, and makes some recommendations for the initial prototypes of the Knowledge Base. The Knowledge Base requirements are examined and then used as criteria for evaluation of the database management options. A mock-up of the data expected in the Knowledge Base is used as a basis for examining how four different database technologies deal with the problems of storing and retrieving the data. Based on these requirement and the results of the evaluation, the recommendation is that the Illustra database be considered for the initial prototype of the Knowledge Base. Illustra offers a unique blend of performance, flexibility, and features that will aid in the implementation of the prototype. At the same time, Illustra provides a high level of compatibility with the hardware and software environments present at the US NDC (National Data Center) and the PIDC (Prototype International Data Center)

ANGULAR CHLOROPHYLL INDICES ESTIMATES DERIVED FROM GROUND-BASED DIURNAL COURSE DATA AND MULTIANGULAR CHRIS-PROBA DATA: TWO CASE STUDIES

ABSTRACT At leaf and plant level chlorophyll indices have shown strong correlations with chlorophyll content and photosynthesis-related processes. However, at canopy level additional abiotic and biotic factors confound the fidelity of these indices. For instance, the Photochemical Reflectance Index (PRI) is known to be sensitive to viewing angles and canopy structure. In this paper we present case studies of two natural canopies at different scales where the influence of sun-target-sensor geometry and canopy structure is inter-compared for a range of chlorophyll indices. In the first case study, surface reflectance was measured in a montane grassland ecosystem located at the Bily Kriz experimental study site (Czech republic) using a stationary mounted AISA (Airborne Imaging Spectrometer for Applications) spectrometer. The experimental set-up resulted in a ground pixel resolution of ~2mm. The effects of changing sun angles on the indices were assessed through a diurnal sampling between 9:00 and 15:00 hrs (local time). Classes of shaded and illuminated photosynthetic (PV) and non-photosynthetic vegetation (NPV) were distinguished per image using a pixel wise classification. The relative contributions of confounding factors as well as the influence of the diurnal variability on performance of the selected chlorophyll indices were evaluated. In the second case study, surface reflectance was measured over an Alpine coniferous ecosystem in the Swiss National Park (Switzerland) using multiangular hyperspectral CHRIS-PROBA (Compact High Resolution Imaging Spectrometer onboard the Project for On-board Autonomy) satellite system with a ground pixel resolution of 17 m. The angular signature of PRI and the structure invariant pigment index (SIPI) was assessed using CHRIS data. Besides, we evaluated the influence of varying tree crown composition and varying viewing angles to the chlorophyll indices with a radiative transfer model FLIGHT. In both cases, the PRI and the green NDVI (gNDVI) responded extremely sensitively to the considered confounding factors at canopy level. The Transformed Chlorophyll Absorption in Reflectance Index normalized by the Optimized Soil-Adjusted Vegetation Index (TCARI/OSAVI), designed to be insensitive to background and LAI variations, responded more sensitively than the conventional NDVI. No certain sensitivity was found for SIPI. The pronounced sensitivity of e.g. PRI and gNDVI, on one hand, and the inconsistency between the chlorophyll indices, on the other hand, erodes the fidelity to use these spectral indices as an effective non-destructive chlorophyll detector

Surface Wave Simulation and Processing with MatSeis

In order to exploit the information on surface wave propagation that is stored in large seismic event datasets, Sandia and Lawrence Livermore National Laboratories have developed a MatSeis interface for performing phase-matched filtering of Rayleigh arrivals. MatSeis is a Matlab-based seismic processing toolkit which provides graphical tools for analyzing seismic data from a network of stations. Tools are available for spectral and polarization measurements, as well as beam forming and f-k analysis with array data, to name just a few. Additionally, one has full access to the Matlab environment and any functions available there. Previously the authors reported the development of new MatSeis tools for calculating regional discrimination measurements. The first of these performs Lg coda analysis as developed by Mayeda and coworkers at Lawrence Livermore National Laboratory. A second tool measures regional phase amplitude ratios for an event and compares the results to ratios from known earthquakes and explosions. Release 1.5 of MatSeis includes the new interface for the analysis of surface wave arrivals. This effort involves the use of regionalized dispersion models from a repository of surface wave data and the construction of phase-matched filters to improve surface wave identification, detection, and magnitude calculation. The tool works as follows. First, a ray is traced from source to receiver through a user-defined grid containing different group velocity versus period values to determine the composite group velocity curve for the path. This curve is shown along with the upper and lower group velocity bounds for reference. Next, the curve is used to create a phase-matched filter, apply the filter, and show the resultant waveform. The application of the filter allows obscured Rayleigh arrivals to be more easily identified. Finally, after screening information outside the range of the phase-matched filter, an inverse version of the filter is applied to obtain a cleaned raw waveform which can be used for amplitude measurements. Because all the MatSeis tools have been written as Matlab functions, they can be easily modified to experiment with different processing details. The performance of the propagation models can be evaluated using any event available in the repository of surface wave events

Imaging the Schwarzschild-radius-scale Structure of M87 with the Event Horizon Telescope Using Sparse Modeling

We propose a new imaging technique for radio and optical/infrared interferometry. The proposed technique reconstructs the image from the visibility amplitude and closure phase, which are standard data products of short-millimeter very long baseline interferometers such as the Event Horizon Telescope (EHT) and optical/infrared interferometers, by utilizing two regularization functions: the ℓ_1-norm and total variation (TV) of the brightness distribution. In the proposed method, optimal regularization parameters, which represent the sparseness and effective spatial resolution of the image, are derived from data themselves using cross-validation (CV). As an application of this technique, we present simulated observations of M87 with the EHT based on four physically motivated models. We confirm that ℓ_1 + TV regularization can achieve an optimal resolution of ~20%–30% of the diffraction limit λ/D_(max), which is the nominal spatial resolution of a radio interferometer. With the proposed technique, the EHT can robustly and reasonably achieve super-resolution sufficient to clearly resolve the black hole shadow. These results make it promising for the EHT to provide an unprecedented view of the event-horizon-scale structure in the vicinity of the supermassive black hole in M87 and also the Galactic center Sgr A*